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ABSTRACT 

We study lensing magnification of source galaxies by intervening galaxy groups and 
clusters using a halo model. Halos are modeled with truncated NFW profiles with cllip- 
ticity added to their lensing potential and propagated to observable lensing statistics. 
We present the formalism to calculate observable effects due to a distribution of halos 
of different masses at different redshifts along the line of sight. We calculate the effects 
of magnification on the number counts of high-redshift galaxies. Using BLAST survey 
data for submillimeter galaxies (SMGs), we find that magnification affects the steep, 
high flux part of the counts by about 60%. The effect becomes much stronger if the 
intrinsic distribution is significantly steeper than observed. We also consider the effect 
of this high-redshift galaxy population on contaminating the Sunyaev-Zeldovich (SZ) 
signal of massive clusters using the halo model approach. We find that for the major- 
ity of clusters expected to be detected with ongoing SZ surveys, there is significant 
contamination from the Poisson noise due to background SMGs. This contribution 
can be comparable to the SZ increment for typical clusters and can also contaminate 
the SZ decrement of low mass clusters. Thus SZ observations, especially for the in- 
crement part of the SZ spectrum, need to include careful modeling of this irreducible 
contamination for mass estimation. Lensing further enhances the contamination, es- 
pecially close to the cores of massive clusters and for very disturbed clusters with large 
magnification cross-section. 
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1 INTRODUCTION 



Gravitational lensing, the deflection of li ght rays from background galaxies by intervening objects, has a number of cos- 
mological and astrophysical applications (Schneider et al.l 19921: iBlandford fe Naravanl 1992 : Naravan fe Bartelmann 19961 : 
Bartelmann fe Schneiderl 1200 ll : ISchneider et al.l 120061 : iHoekstra fe Jainl l2008h . In the weak lensing regime, small distortions 



of source galaxies can statistically constrain cosmological parameters related to structure formation and theories of gravity. 
In the strong lensing regime, where light is deflected by rare massive galaxy clusters, it is also possible to infer detailed 
information on cluster profiles, though with less statistics. Although lensing conserves surface brightness of sources, it changes 
their observed fluxes and sizes; e.g. background galaxies have their fluxes magnified in the line-of-sight of massiv e halos cores. 
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rovid e d one can detect them in pure samples, and measure their masses and redshifts precisely 



20071 : IRozo et al.ll2007l , l2009h . Techniques for cluster detection and mass measurement include 



the counting of optical galaxies, the measurement of their lensi ng signal, their X-ray temperature or flux and their Sunyaev- 
Zel'dovich flux decrem ent /increment (SZ effect). The SZ effect I Sunvaev fc Zeldovich|[l972l : lltoh et al. 1998 : Birkinshaw 19991 : 



Carlstrom et al.ll2002f ) results from the up-scatter of cosmic microwave background (CMB) photons by the hot electrons in the 
intra-cluster medium, which shifts the underlying spectrum and causes a decrement/increment at wavelengths longer/shorter 
than ~ 1 mm. This allows for the detection of clusters independently of their redshifts and the measurement of their masses, 



* mlimaOsas.upenn.edu 
© 0000 RAS 



2 Lima, Jain and Devlin 



as long as the effects of possible contaminants, such as radio sources and far-IR/submillimeter galaxies can be removed or at 
least well understood. 

The submillimeter galaxy population (SMGs) consists of dusty high-redshif t galaxies with high star formation rates 
( Blain fc Longairl 1 19931 ; fearger et all Il998l ; iHughes et all Il99sl ; iBlain et all I2OO2I ). The energy output of the star forming 
regions heats up dust grains, which emit a modified blackbody spectrum that peaks around 0.1 mm in the rest-frame. SMGs 
can produce a non- negligible point-source contribution to CMB anisotropics ( Scott fc White! 1999h . When positioned on the 
line-of-sight of cluster cores, these galaxies can contaminate the cluster SZ flux, fillin g up the decrement and enhancing the 
incre ment, resulting in mis-estimation of the cluster mass from SZ measurement s ( Knox et al. 20041; White fc Maiumdar 
2004). In addition, galaxy clusters can mag nify the fluxes of the SMGs (jBrainll 19961 . Il997l . ll999l ; iPaciga et al.ll2008r i. changing 



their counts distribution and further enhancing the contamination of the SZ signal. 

In this paper, we study the effect of lensing magnification of SMGs by foreground clusters from a twofold perspective. 
On the one hand, galaxy clusters magnify the fluxes of SMGs, changing their intrinsic counts distribution, but potentially 
facilitating their detection and follow-up studies. On the other hand, the magnified fluxes of SMGs can cont aminate unresolved 
SZ clu sters and may need to be understood in order to provide clean SZ mass estimates. For instance, iDiego fc Partridge 
( 20091 ) find that including point sources in the modelling of cluster masses results in better agreement between x-ray and SZ 
mass estimates. Follow-up observations of unresolved clusters provide a means to tackle these issues, allowing for separation 
and removal of SZ contaminants and detailed studies of their properties. 

We start in § [2] describing the properties of the dark matter halos used to represent the lens galaxy clusters. In §[3] we 
briefly review the basic formalism to analytically compute lensing magnification by massive halos and in § [4] we describe the 
lensing probability and some of its applications. Our main results are shown in § [5] where we estimate the effects of lensing on 
the SMGs distributions (§ 15. 1|) and show the contamination effect of the (lensed) SMG population on the SZ flux of clusters 
(§E21). We conclude in §H 

Throughout, we assume a fiducial cosmology for a flat univer se with parameter v alues based on the results of the Wilkinson 
Microwave Anisotropy Probe first year data release (WMAP1) jSpergel et all \200j ) The cosmological parameters (and their 
values) are the normalization of the initial curvature spectrum <5<;(= 5.07 x 10 -5 ) at k = 0.05 Mpc -1 (corresponding to 
ag = 0.91), its tilt n(= 1), the baryon density relative to critical Qbh 2 (= 0.024), the matter density f2 m /i 2 (= 0.14), and two 
dark energy parameters: its density J7de(= 0.73) and equation of state w(— —1), which we assume to be constant. 



2 HALO PROPERTIES 

In this section we present the halo properties and scaling relations as sumed throughout to estimate lensing effects of galaxy 
clusters. We take the NFW profile prescription for dark matter halos (jNavarro et al.l ll997'l 

ps 



p{j ' (cr/r vir )(l + cr/r vir ) 2 ' 
where r v j r is the virial radius and the halo concentration c is given by a fit to simulations from Bullock et al. 1 200 ll ) 

9 (M vi -\- 013 

c(M vil ,z) = 



1 + z V M. 



(1) 



(2) 



with M* such that cr(M») = S c . Here a 2 (M) is the variance of the linear density field, defined in Eq. (O below, and 8 C is 
the linearly extrapolated density contrast threshold in spherical collapse. Even though S c has a small redshift and cosmology 
dependency, here we take it to be fixed at its value in a fi m = 1 universe, i.e. S c — 1.686. The halo virial mass M v i r is 



Mvir 



drAirr p(r) 



r 



(3) 



where / = (ln(l + c) — c/(l + c)) . The virial radius r v i r can be computed from the virial overdensity A c relative to critical 
density 

A3A/ V ir / . \ 

c = 1--. 3-, (4) 

where p CT it(z) = p C rit,oE 2 (z) is the critical density at redshift z, p cr it,o is its value at z = and the scaled Hubble parameter is 



We take A c from a fit to simulations of iBrvan fc Norman jl998h for flat cosmologii 



(5) 



A c = 18tt 2 + 82x - 39a: 2 



(6) 
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Figure 1. (Left): Convergence k, shear 7 and magnification (iasa function of angular separation 6 from the halo center. Here the lens 
mass and redshift are M v ; r = 10 14 h~ 1 Mq and = 0.2, the source is at redshift z s = 1.0 and the halo is assumed to have a NFW profile. 
(Flight): The virial radius 8 vlr in angular units as a function of Z\ for different halo masses. 



where x = u m (z) — 1 and iOm(z) — fi m (l + z) i /E 2 (z). One can similarly define the overdensity relative to the mean matter 
densit y A v i r = A c /w m (z)- Finally for the halo distribution in mass and redshift we take the description of ISheth fc Tormen 
( 19991 ) 'or the comoving differential number density of halos per logarithmic mass interval 

dn p m dv 



dlnM vir M vir dlnM v i r ' 
where v = 5 c /cr(M v - a ) and 



uf(v) = Ay ^av 2 [l + {av)~ v ] exp[-av /2] . (8) 

Here a 2 (M) is the variance of the linear density field in a top hat of radius r that encloses M = 4-7rr 3 p m /3 at the background 
density 

« 2 (r) = J ^\W(kr)\ 2 P^k), (9) 

where Ph{k) is the linear power spectrum and W is the Fourier transform of the top hat window. The normalization constant 
A is such that f dvf(y) = 1 and we take the parameter values p = 0.3, a = 0.75. 



3 LENSING MAGNIFICATION BY MASSIVE HALOS 
3.1 Axially Symmetric Lenses 

Given the 3d density profile p(r), we can compute all lensing quantities of interest in the lens plane by performing successive 
numerical integrations and differentiations. For the NFW profile, most quantities can be computed analytically, speeding up 
calculations that use the lensing properties. The projected density field E is obtained integrating over the parallel coordinate 
r« = x m the position vector decomposed as x = (ru,rj_), where the perpendicular coordinate r± — Da(x)^i Da is the 
angular diameter distance and 6 is the angular coordinate in th e lens plane. Here and throughout, all distances are comoving. 
For the NFW profile the projected surface density is given by (jTakada fc Jainll2003ah 



e(<?) = p" dr„ P (r) = ^^-Hce/e^) , (10) 
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Figure 2. (Top Left): Magnification (iasa function of angular separation 6 from the cluster center, assuming a NFW profile, for different 
cluster masses with z; = 0.2 and z B = 1.0. (Top Right): Same but for z; = 0.5 and z s = 2.0. (Bottom Left): fi evaluated at 6 = 0.5 arcmin 
as a function of lens redshift z; for z s = 1.0 at different cluster masses and for z s = 2.0 for cluster mass M v ; r = 10 14,2 /i~ 1 Mq. 
(Bottom Right): fi(8 = 0.5arcmin) as a function of source redshift z s for zj = 0.2 at different cluster masses (upper panel) and for 
M v ; r = 10 14 2 h~ 1 Mq at different values of zj (lower panel). 



where 



F(x) 



\J c2 ~ x2 j i 

(l-x 2 )(l + c) + (1-x 2 ) 3 / 
3(l + c ) [ L ^ c+lj ' 



■arccosh 



(* = !) 



(11) 



' (l-z^l + c) ( I 2_ 1) 3/2 ' 



0, 



(a; > c) 



(12) 



and 6'vir(x) — / Da{x)- The convergence field n is defined in terms of the critical density E cr i t by 

( > ~ E cr5t ' ^ ctit - A-kGD a { X )D a (xs-x)' 
where a = (1 + is the scale factor. The convergence measures isotropic light distortions and is related to the projected 
lensing potential tp via the Poisson equation \7 2 ip = 2k. Anisotropic distortions are measured by the complex shear field 
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7 = 71 + 172, whose amplitude for an axially symmetric mass distribution is simply related to the convergence by 



7 (e) = r(<0)-«(0), 

where the average convergence field up to 9 is 

r.0 



_L_ 



d9'2%9'K,{9) . 



(13) 



(14) 



7(*) = 
where 



Like the convergence, the shear can be computed analytically for the NFW profile and is given by (jTakada fc Ja in 2003b) 
_ M vil fc 2 G(e0/0 vir ) 



(15) 



G{x) = < 



x' 2 (l + c) 

1 

3(l + c) 
1 

f- 1 



(2-x 2 )^/c 2 -x 2 



(llc+10) yV 2 
1+c 



(2-x 2 )yJ c 2 -x 2 




(16) 



Finally, the magnification /i measures the total flux amplification as well as the increase in angular size of source galaxies, 
and is given by the inverse determinant of the Jacobian transformation between source and image angular coordinates. It is 
expressed in terms of the convergence and shear as 



H(0) = 



I -«(*)]' -| 7 (0)|a 



(17) 



In Fig. [T] we show various lensing observables, assuming a spherical NFW profile lens halo of mass M vir = 1O 14 /i _1 M 
at redshift zi = 0.2 and a source galaxy at z 3 — 1.0. The convergence tends to dominate over the shear in the halo core 
whereas the shear dominates in the outer parts. Likewise the magnification is close to 1 (no magnification) far from the halo 
and rises to significant values close to the cluster center. Also shown is the virial radius in angular units 9 va as a function of 
halo redshift for different halo masses. Halos of interest have typical angular sizes of the order of few arcminutes. 

In Fig. [21 we show the magnification fi as a function of angular separation 9, lens redshift Zi and source redshift z 3 
for different cluster masses. As the halo mass increases the magnification is significantly enhanced. The two spikes in each 
magnification curve at the top panels of Fig. [2] represent critical curves (tangential and radial) where the magnification is 
formally infinite and correspond to the two solutions of the quadratic equation |1 — k\ 2 — ■y 2 . 

Magnification increases with lens redshift, reaching a maximum, and then decreases as the lens approaches the source; 
that reflects the lensing efficiency of the critical surface density terms. Magnification also increases with source redshift, until it 
reaches a plateau; therefore one expects roughly similar lensing magnification effects from sources at sufficiently high redshifts. 

In what follows, we use the letters k, 7 and fi to denote full lensing quantities, which incorporate halo ellipticity (§ 13.21) . 
The spherically symmetric results of this section will be denoted by the same letters with a tilde (e.g. k). 



3.2 Elliptical Lenses 

Ellipticity can be introduced in a number of ways into halo profiles. The most obvious choice is to introduce it directly into 
the 3d density profile p(r) and project it to obtain E. Another choice is to introduce ellipticity in E directly. The computation 
of the lensing quantities can then be obtained numerically by performing integrations (of E) to obtain deflection angles, and 
differentiations (of deflection angles) to obtain shear and magnification fields. In special cases of known potential-density pairs, 
these quantities can be computed analytically. 

Another possibility is to introduce ellipticity in the two-dimensional projected lensing potential <p, whose second derivatives 
give the lensing quantities more directly. The latter approach may produce density profiles with dumbbell shapes for high 
values of ellipticity , but produces phy sical magnification fields for low/intermediate ellipticities close to those expected for the 



majority of halos (jjing fc Sutdl2002r ). Moreover, this approach leads to simple analytical expressions that can be expressed 



in terms of the usual results of the spherically symmetrical case and, for this reason, this is the choice we implement. 
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Figure 3. Maps of magnification fj, for different values of ellipticity with z\ = 0.5 and z B = 2.0 and halo mass M v j r = 10 15 /i _1 Mq. Here 
e = 1 — b/a, where a and b are the projected major and minor ellipse axes in the lens plane. (Top panels): Ellipticity e v is introduced 
directly in the potential profile using the analytical expressions derived in this paper. (Bottom panels): Ellipticity e K is introduced in 
the density profile and propagated to the deflection angle and magnification numerically using the WSLAP code of Diego et al (2007). 
Notice that e v = 0.2 produces roughly the same magnification map as e K = 0.4. The field of view is 0.6 arcmin in size. 



3.2.1 Elliptical Potential 



Given an axially symmetric lensing potential ip(0), we follow iMeneehetti et all |200ot ) and obtain the elliptical generalization 



with major axis along the 62 direction by substituting 6 



'I by 



+ 6i{l-e v ), 



(18) 



b/a with a and b being the major and minor ellipse axes. Our approach is similar to that of lGolse fc Kneib 
1 20021 ). 1 hough our NFW profiles are truncated at the virial radius as opposed to extending to infinity. We can then compute 
the deflection angle components, their derivatives and the lensing quantities of interest for elliptical lenses. In Appendix A, 
we express the results for re, 71 and 72 in terms of the spherically symmetric case, with the effects of ellipticity described 
analytically. Using the expressions given in Eqs. lAlfl IA12l and lA13l we can compute the magnification defined by Eq. (|17[) . 

In Fig. [3] we show magnification maps in the 9\ x 62 plane, for different ellipticities. The top row shows maps derived using 
our model in which ellip ticity is introduced in the lensing potential. The bottom row shows maps derived using the WSLAP 
code ([Diego et al.ll2007l ). where ellipticity e K is introduced in the 3d density profile, which is then projected. Notice that both 
agree when e v = e K = and that the ellipticity introduced in the potential causes larger changes in the magnification map. 
This is expected, since the second derivative of the potential is related to the density - the quantitative connection between 
the two is presented below. 

In Fig. [4] we show the magnification /x as a function of angular separation along the two directions 6\ and 62 for different 
values of e v . In the 62 direction of the major axis, the tangential caustic moves to larger angles whereas the radial caustic 
moves to smaller angles. In the 9i direction much less dramatic changes happen, as also seen in Fig. [3] 
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Figure 4. (Left): Magnification fi as a function of angular separation from the cluster center in the 0i-dircction (82 = 0). (Right): Same 
but in the 02-dircction (9± = 0). We assume an elliptical NFW profile, with different values of ellipticity, Zi = 0.5 and z B = 2.0. Here 
e v = 1 — b/a, where a and b are the major and minor ellipse axes and ellipticity is introduced in the lensing potential in the 02-direction 
with the prescription of Meneghetti et al. (2003). 



3.2.2 Elliptical Density versus Elliptical Potential 

Simulations typically predict ellipticities in the halo density profile. Therefore it is interesting to relate the ellipticity introduced 
in the potential as in the previous section, denoted e v , to the corresponding effective ellipticity in the surface density profile, 
denoted e K , introduced directly in the convergence map instead of the potential, i.e. by substituting 8 2 — 8\ + 8\ in k by 



(19) 



(1 

In Fig. [5] we show convergence maps for different ellipticities, similarly to Fig. [3] Again, ellipticity e v is introduced in the 
lensing potential (top row) and e re in the density (bottom row). Fig. [6] shows one way to relate the two kinds of ellipticity. 
Enforcing the match ^(81,82) — «(#*), it shows the relationship between e« and e v . Here ^(81,82) is the convergence that 
results from adding in the lensing potential and k is the spherically symmetric convergence. Specifically we match the 
value of the convergences at the arbitrary matching point (81 — 0,82) for various choices of cluster mass, lens and source 
redshift and matching coordinate 82- The relation is quite insensitive to these parameters and also to the value of 82 chosen 
for the match. Overall the relation e v — 0. 48e K seems to hold rel atively well for e K < 0.5 in the range of parameters considered 
here. This is consistent with the results of Golse fc Kneibl 1 20021 ) for the extended NFW profile, where they find e K ~ 2e v for 



e v < 0.25. Our results can be used to relate the typical values of e K from simulations to the more convenient values of e v 
used to analytically model the effects of ellipticity. 



4 LENSING STATISTICS 
4.1 Halo Counts 

Given the differential comoving number density of lens halos rfn/rfln A/ V ir, one can estimate the cluster number counts dN in 
the redshift/mass range dzidln M V - 1T and solid angle Afl in the sky as 

N t 7-1 dV(zi) dn(zi,M v ir) 

dN(zi,M vir ) = dzidlnM viv , v \ ' — '- , (20 

dzi dmAivir 
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where dV/dzi — AflD\/H is the comoving volume element, Da(zi) is the comoving angular diameter distance and H(zi) is 
the Hubble parameter. The differential halo density n and counts -/V above a mass threshold M t h are given by 

, x r w dn(zi,M vil ) 

n(2i) = J Mtb dlnMvir dlnM vir ' (21) 

— ; = Ai2 E77 ^ n \ Zl > ■ 

dzi H(zi) 

4.2 Lensing Cross-Section 

Each halo produces a elliptical "ring" region around it, with area AQ^(zi, z„, M, e v , fimin), such that inside this area the 
magnification is larger than fi m in- This area is an effective cross-section for lensing statistics. Note Af2 M is the area with 
/i > l-imin in the source plane 

o2 f d8 2 



AQ ll ( t x min )= / dF= -fx, (22) 

where we used the fact that the magnification is precisely the Jacobian of the transformation between image and source 
coordinates /i = d8 2 /df3 2 . A reference value for Af2 M is the projected cluster surface area defined by its virial radius in the 
image plane Afi v i r = n8 2 iT . For large values of /i we expect the magnification region to be well inside the halo core and 
therefore Af2 M <g Aflvir, whereas for /x — > 1, Af2 M becomes formally infinite. For a given halo, we estimate Af2 M by evaluating 
H in a square grid of 100 x 100 points with side 29 v \r- We then multiply the fraction of grid points with fi > fi m in by the area 
4#vi r ; in this fraction, each point is weighted by l//j,(0) at the grid. 

Note that the way we estimate Ar2 M is only strictly accurate if Af2 M < 48 2 iT . This is not true for sufficiently low fi, where 
A£7 M formally covers the whole sky. In principle this estimation can be improved by adaptively surveying a radius larger than 
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Virial and lensing cross-sections (in rad 2 ) as a function of halo mass, redshift, ellipticity and minimum magnification. (Left): 
and An v ; r = 7r0 2 ir as a function of halo mass M v ; r for z\ = 0.5, z s = 2.0, e v = 0.2 and fimi n = 10. (Middle): Same but as a 
2; for M v ; r = 10 15 /i _1 Mq. (Right): Af2 M versus e v for fi m i„ = 10 (upper panel) and versus /i for e v = 0.2 (lower panel). 



2# v i r around the halo. For the goal of estimating the lensing probability, this is not necessary as discussed in the next two 
sections. 

In Fig. [7|we show Afi M and Af! v i, as a function of halo mass, redshift, ellipticity and minimum magnification. Since 
r v i r ~ M V 3 , we have that Af2 v i r ~ . For this high value of /U m i n , AS7 M is many orders of magnitude smaller but growing 
faster with M vir . Since r vir ~ [A c (z)E 2 (z)}' 1/3 , we have that Af2 vir ~ A c (z) _2/3 £(z) _4/3 Da(z)~ 2 . The redshift dependence 
of Afivir is dominated by that of Da, since A c (z) and E(z) change less. In particular, it increases rapidly at low redshifts, 
where Da ~ > 0. On the other hand Afi M , which shows a similar trend at intermediate redshifts, goes to zero as zi — > or z s , 
where the lensing efficiency vanishes. 

As also apparent from the magnification maps, there is a relatively weak dependence on e v , where ellipticity mainly 
distorts the critical curves, but do not change their area. A strong dependence on fj, m in can be seen; the largest areas come 
from Umin ~ 2 even for the most massive halos. We will see below that once we sum over halo mass, the contribution of high 
magnification (/i £ 2) regions to observable number counts is very small. 
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4.3 Sky fraction with large magnifications: Lensing optical depth 

We can estimate the fraction of the sky with fi > fimin due to all halos above a certain mass and redshift range as 

ffj, = dzi^^- r rflnAfvir r dz s P{z s ) f\e v P(e v )A^(z l ,z s ,M vir ,e v ,^ mm ) dn j^\^ v[l \ (23) 

Jo H ( z l) JMa, Jo Jo dlnMvir 

where P(z s ) is the source redshift distribution and -P(e^) is the distribution of halo ellipticities. This quantity is the opti- 
cal depth for lensing magnification and determines the probability that a given source galaxy /population is magnified by 
intervening lens halos. This definition of / M accounts for overlapping angular regions, which are counted multiple times in 
the integrals in mass and redshift. However it neglects multiple lens events, such as multiple images, since only one event is 
counted per halo. Therefore, strictly speaking, / M is the sky area with large magnifications only for / M <C 1, i.e. fj, ^> 1, where 
no such overlappings occur. 

Note also that / M can be, and in fact is, larger than 1 at regions with ji ~ 1, i.e. such regions cover the sky many times. 
In fact, even AS7 M itself for a given cluster can be as large as the whole sky for fj, ~ 1. How large becomes as ^ approaches 
1 depends on the top value we allow AS7 M to have, which in our case was set as 4# 2 . ir and causes a saturation in / M at ~ 27 
for fj, ^ 1. Even though this top cut-off in Afl M affects our computation of for very low values of fx, it does not significantly 
change the quantity we are really interested in, the lensing probability P(> ^i) described in the next section, which in any 
case is close to unity for sufficiently large / M . 

We can also define a fraction / v i r , similarly to / M , replacing AS7 M — > Af2 v i r . It turns out that the area of halos defined 
in this way can be larger than the sky area (/ V j r > 1) if one includes halos of sufficiently low mass. However, those halos are 
unable to produce large magnification areas and therefore / M remains smaller than 1 for sufficiently large values of Umin- 

For fixed source redshift and cluster ellipticity (in which case P(z s ) and P{e v ) are delta functions) we have 

du _ r dlnMra%(ziiMviriMmm) ^ii, (24) 



dUr _ r dlnMvirA ^ r( , liMvir) *^M. (25) 



d Zl H{ Zl ) J Mtb - „ dlnMv 

In Fig. [5] we show the sky fraction dfvir/dzi and df^/dzi as a function of zi for fixed z s = 2.0, e v = 0.2 and fi m in = 10. 
The latter fraction does not change much for M t h < W 12 h -1 M@ since these halos, though numerous, are unable to produce 
large magnifications. Very massive halos produce large magnifications but their abundance is exponentially suppressed. 

In the left panel of Fig. [9]we show the sky fraction as a function of fimin for Mth = 10 12 M@, z s = 2.0 (left). The 
fraction increases sharply as \i m in approaches 1. The curves with different values of e v show that for large magnifications, / M 
can be significantly higher for large ellipticities. 
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Figure 9. {Left): Sky fraction / M with magnification it > fJ, m in as a function of (J, m in Ior halos with masses M V \ T > 10 12 /i — 1 Mq and 
different ellipticities e v . (Right): Lensing probability P(fJ.) to have magnification along the line of sight up to z B = 2.0 from all intervening 
halos, assumed to have e v = 0.2 (dotted line) and = 0.4 (dashed li ne). Also shown is the lensing probability obtained from N-body 
simulations (thin solid line) for sources at z s = 2.1 l|HiIbert et al.|[2007f ). 



4.4 Lensing Probability 

Given the lensing optical depth / M we obtain the lensing probability as 
P(> n) = 1 - e~ /f ' . 



(26) 



Note that P(> (J,) ~ / M for / M < 1 or large magnifications, i.e. small optical depths are themselves the lensing probability. On 
the other hand, P(> /i) — ► 1 for / M — > oo or small magnifications. Our top cut-off value of 48^ iT for Af2 M implied a saturation 
of / M ~ 27 for ii ~ 1. At this value of / M , however, P(> fi) is already very close to unit for all practical purposes and therefore 
is not affected significantly by the cut-off. 

Note P(> n) is the probability that a source galaxy is magnified by more than fi and therefore it can be expressed in 
terms of the probability density P((J,) as 



P(>fi) 



P(fi)dfi. 



(27) 



Alternatively P(fJ,) can be obtained by P(fi) = —dP(> fi)/dfi. This statistics can then be incorporated in various applications 
to estimate the effects of lensing magnification. An interesting property of the lensing probability in the source plane is that 
P(> At) oc and th erefore P(fi) oc 1/ii 3 for fi 3> 1, as can be shown in particular cases and argued to be true in general 

( Schneider et al.lll992l ). 

In the right panel of Fi g. [9] we show Pin) for z B = 2.0 and e v = 0.2, 0.4. We also show the magnification probability 
for sources at z s — 2.1 from iHilbert et al.l l|2007l ). obtained by ray-tracing high- resolution N-body simulations. Notice that, 



for /j > 1, P(> A*) ~ /m k A* 2 as shown in the left panel of Fig. [9] and it follows that P(ai) oc fj, 3 as expected. Our 
semi-analytical estim ate of P(p,) with e„ = 0.4 compares relatively well with the probability obtained from ray-tracing on 



N-Body simulations (jHilbert et al.ll2007l : iPaciga et al.ll2008l ) at large magnifications. 

Even though e v = 0.2 is in better agreement with the average values of halo ellipticities Uing fc Suto 2002), the fact 
that our model does not include halo substructure and other realistic effects such as other halos along the line of sight 
that enhance magnification appears to be incorporated by artificially increasing the halo ellipticity to e v = 0.4. In addition, 
the semi-analytical model underestimates P(/i) compared to the simulation results at low magnifications (p, < 1.0). That 
is expected since our model only includes a one-halo term and does not properly account for the effects of voids, where 
highly de-magnified regions sit. That feature is well captured by simulati ons and can potentiall y be obtained in an improved 
version of our model th at includes compensated halos, a two-halo term (|Coorav fc Shetbll200 2) as well as halo substructure 
i Meneghetti et all 20051). In our estimations of lensing results on high-redshift galaxy populations, we use P(/i) from N-body 
simulations (|Hilbert et al.ll2007l ). 
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Finally, note that P((J.) denned above averages over clusters in the line of sight up to the source redshift z 3 . When we 
are interested in the effects within a specific cluster solid angle, we must use only the region around the cluster to define the 
lensing probability. For instance, if we want to estimate lensing effects within the virial radius of a specific cluster, we use a 
probability defined simply as P(> /z) = Afi M /Afi v ir, where these quantities are calculated from our semi-analytical model 
described above. 



5 OBSERVABLE EFFECTS OF LENSING MAGNIFICATION 



5.1 Counts of Submillimeter Galaxies 

Submillimeter galaxies (SMGs) are dusty galaxies, typically optically obscured but visible in wavelengths below 1mm due to 
the thermal spectrum re-emitted by the dust grains after absorbing radiation from the star-forming regions. The simplest 
prescription for the Spectral Energy Distribution (SED) of SMGs is it is given by a blackbody spectrum £>„ at a high 
temperature T ga i, modified by an emissivity e tx 



eB v 



,3+13 



exp 



hv 



- 1 



.3+13 



exp 



(z//GHz) 



21(T ga l/K) 



- 1 



(28) 



Apart from an overall normalization factor, the 2 free parameters (5 and T ga i as well as the galaxy redshift z 3 must be 
known. Typical fits to observed galaxies give high redshifts of z 3 ~ 1 — 6, values of ft in the range 0.5 — 2 and temperatures 
of ~ 15 — 40K. Specifically, in our fiducial prescription we assume that all galaxies are at z 3 = 2 (though we also show 
some results for z 3 — 3.0 and 4.0), and take T ga i = 30K and ft — 0.7. In the Rayleigh- Jeans regime of low frequencies, the 
temperatur e is un i mportant and S,; ED o c u 2+lS . The value of /3 = 0.7 is consistent with the spectral index a — 2 + /3 used by 



Knox et all ((200J); IWhite fc Maiumdarl ([2004]). Using the SED prescription one can scale fluxes at one observed frequency v 
to another frequency of interest v' ', e.g . in the SZ flux range . As we will see, ft — 0.7 properly accounts for extrapolated fluxes 



in the empirically motivated model of lLagache et al.l (|2004Q for the SED of SMGs at different frequencies. 

The count dist ribution of SMGs have b een measured at va rious wavelengths, including at 250, 350 and 500/^ m by BLAST 



Devlin et al 



Pereraetal 



20091) . at 850>m by SCUBA JCoppin et alfood). at 1 100/i m by BOLOCAM (lLaurent et alJbood ) and AzTEC 



20081 ) and at 1200/itm by MAMBO (jGreve et al 



20041 ). See lPearson fc Khar] (|200g| ) for a recent compilation of 



these measurements, a model for fitting them and prospects for future measurements. 

In our standard prescription we take the differential counts dn/dS v of the submillimeter galaxy population as measured 
by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) at A = 500/^m (y ~ 600GHz), and scale various 
quantities of interest to other frequencies according to the SED prescription. When studying the effects of lensing on the 
properties of SMGs, we also consider the example of a Schechter model, in which the exponential steepness of the counts make 
the lensing effects more dramatic. 



5.1.1 Magnification and number counts 

As a result of photon and energy conservation, the surface brightness of galaxy sources, defined as the flux per unit solid angle, 
is conserved by gravitational lensing. Since magnification, by definition, increases the solid angle Q of sources by a factor /it, 
it has to also increase their intrinsic flux S, effectively lowering the observed threshold and increasing the number of sources 
available 

S ->■ S ohB = vS, (29) 
dii -> dQ obs = fidSl . (30) 

As a result, the intrinsic number density of a source population is modified by lensing m agnification. For a given magni- 



fication fi, the intrinsic differential number density dn/dS is modified as (see e.g. Eq. (11) in lRefregier fc Loebl (|l997r i) 
dn 1 dn ( SobA 

ds ~* M 2 ds { (i J ■ {6 > 

The 1/fj, 2 factor comes from transforming the angle differential dfi (implicit in the above equation) and the flux differential 
dS into their observed counterparts. The additional change is from the fact that the observed flux S b s corresponds to a true 
flux of S = S bs/p- The intrinsic differential number density can be interpreted as the probability density P(S) of having 
intrinsic flux S from a random source. Likewise, the observed modified distribution represents the probability density P(S bs\li>) 
of having an observed flux Sobs given a magnification /i. In order to obtain the probability P(S bs) of having Sobs irrespective 
of magnification, we must multiply by P(fi) and integrate over all values of \i. 

Given the differential number density, a number of quantities are of interest, such as the cumulative number density, 



© 0000 RAS, MNRAS 000, 000-000 



Lensing Magnification: Implications for Counts of Submillimeter Galaxies and SZ Clusters 13 



the cumulative number counts N(> S) in a solid angle AQ, the average total flux of the galaxy population (5*) gal and its 
associated Poisson variance (5 2 ) gal 

n(>3) = J^dS', (32) 

N(> S) = Afi n{> S) , (33) 

S = Anj S^dS, (34) 

S 2 = Afi f S 2 ^dS. (35) 
J dS 

Ignoring the clustering of the galaxy population, which adds extra noise to the average flux, the expected fluctuation in 
the average flux of galaxies within solid angle AS! due to Poisson noise in the counts is given by 

a g *\S) = VW. (36) 
The average flux in one frequency v can be scaled to v' given a SED prescription Sf! ED 
r dn S' S , ED r dn S^F - 

s »' = An J s ^ dS " ~ s|^ An J Sv ik dSv = 4^ s - (37) 

and a similar rescaling can be applied to cr ga (Sv). Lensing magnification modifies these various quantities to their observed 
values according to the effect of a given magnification fi and the lensing probability P((J,) of actually having that magnification. 
Since we define P((J,) in the source plane, we may identify angular averages in this plane with averages over -P(/i). For a variable 
X, we have 

^ J dQX = j d»P(fj,)X . (38) 
This implies 

< M > = J duPMfi = ^ . (39) 
Using the relation between averages in the observed (image) and intrinsic (source) plane, i.e. 

Aib / dn ^ x = WKnI » dnx = <J> / d ^ p ^ x ■ ( 4 °) 

and the fact that observations can only average in the image plane, we obtain the observed lensed quantities in terms of their 
unlensed counterparts 

dn ohs (S bs) _ 1 f djl P{^) dn ( ffobs A 



dS obs (n) J fi dS \ (i 

n obs (> Sobs) = A [dfiPQi)n (> — ) , (42) 



iVobs(> Sobs) = TT / dnP(n)N (> , (43) 

Sobs = S, (44) 



S 2 obs = ^S 2 , (45) 



where 

<^> = | M Q P( M ) . (46) 

When considering effects of all halos in the line of sight of source galaxies at large sky patches, we use the P(fi) defined in 
8 )4.41 When considering the effect of the source population on a single cluster, we use P(/j,) defined within the cluster radius. 
Even though an intrinsic solid angle Afiint corresponds to an observed solid angle Af2 b s = (/i)Afii n t, here we assume that 
the angle implicit in Eqs. (|43I45|) is simply fixed by the observation, i.e. Af2 b s = Afi, where AQ was used to define unlensed 
quantities Eqs. (|33I35[) . Notice that Eq. (|44[) reflects the conservation of total surface brightness by lensing and along with 
Eq. (|45[) . is only valid when these quantities are integrated over all fluxes. If we introduce lower or higher limits of integration, 
we need to integrate over P((J,) as for the other quantities. In that case Eq. (|44[1 is no longer true, i.e. the lensed and intrinsic 
brig htness differ. 

Paciga et al. ( 20081 ) considered the lensing effect on the distribution of SMGs detected at 850 ^m. However they did not 



include the factors of l/(/u) and 1/fj, in the expression for the differential counts Eq. (|41[) . That caused an overestimation of 
the magnification effects, which are highly sensitive to the steepness of the differential counts. 
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Figure 10. (Left): Intrinsic differential number density dn/dS from BLAST 500^tm data (points and solid black line) as well as observed 
lensed counts dn obs /dS b s at different source redshifts (dotted and dashed lines). Here the observed BLAST data is assumed to be the 
true underlying distribution at different redshifts, which is then lensed by intervening halos. (Right): Cumulative number density for the 
same unlenscd/lenscd BLAST distribution. 



If the intrinsic differential number density can be parameterized as a power law with index a near the flux threshold, 

^ oc fiT- , (47) 

then the above expressions become simply their intrinsic counterparts rescaled by moments of the magnification probability 
P(ji), e.g. 

dn ohs (S obB ) (/i" -1 ) dn(S ohE ) 

dS ohs fa) dS ' [ ' 

and similarly for n Q b s (> .Sobs) and other observables of interest. Notice that since P(n) oc 1/fJ. 3 as fi — > oo, the observed 
quantities diverge for large enough values of a, unless we impose an upper integration limit defining a maximum magnification 
or a cut-off in the counts distribution. In practice there is always a maximum magnification imposed, for instance, by the size of 
the source galaxies. Moreover, a constant power law index does not hold for counts over all fluxes. For a non-evolving population 
in Euclidean flat space a = 2.5; larger /smaller values of a indicate a luminosity or density that is increasing/decreasing in 
time, i.e. non-trivial evolution of the galaxy population. 



5.1.2 Number Counts from BLAST 



In Fig. [10] we illustrate lensing effects on the distribution of SMGs for the BLAST 500/^m data. We show the measured 
BLAST distributions, which we take to be intrinsic as opposed to observed for illustrative purposes. Shown in Fig. [10] are the 
differential and cumulative number density, as well as the corresponding lensed distributions, assuming sources at different 
redshifts. We spline interpolate the original BLAST points to obtain a smooth distribution at arbitrary values of flux density. 

Obviously, the measured distributions already have lensing effects convolved, which ideally one would try to deconvolve to 
obtain the true distributions. If the lensing effects are relatively small, as in this case, the above assumption is self-consistent 
since the observed distribution is in fact already close to the true one. The data indicates that the observed differential counts 
have a index a ~ 2 until S ~ 10 mjy and a ~ 4 beyond this flux. Notice that for ot > 3 the differential number density 
increases significantly (see Eq. I48p and the relative differ ence between the lens ed and unlensed curves increases at the bright 
end by up to ~ 60%. The empirically motivated model of lLagache et al.l ([2004]) predicts a distribution at 500^im that is quite 
close to the one measured by BLAST, suggesting that the observed slopes are already close to the true ones. Howe ver, in order 



for the total brightness inferred from the BLAST points not to exceed the Cosmic Infrared Background (CIB) (|Dwek et al 
199Sl ; Fixsen et al. 19981 ; Smail et al. 1997 ; Hughes et al. 19981 ). it is necessary to impose a sharp cut-off flux of S cut = 4.6 mjy 
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Figure 11. Intrinsic differential density dn/dS from BLAST 500^tm data (points and thin solid like) as well as observed lenscd counts 
d-n bs/d-S b s at different source redshifts (dotted and dashed lines). Here the true underlying distribution is assumed to be of a Schechter 
type (thick solid line), which is then lensed by intervening halos into the observed BLAST distribution. The underlying distribution fits 
the middle BLAST points and, when lensed, is roughly consistent with the brightest point. 



( Devlin et al 2009) . In § 15.1.31 we consider a case that does not require such cut-off, where fluxes decay at the faint end and 
have much steeper intrinsic slopes at the bright end. 

Here we have integrated up to a maximum magnification of /i max = 100. However, the results remain nearly the same even 
for /i max = 10, reflecting the rarity of high magnification events displa yed in Fig. [9] Therefore, o ur results are insensitive to 
a magnification cut-off introduced by finite size of source galaxies (e. Jxakada Ha m J ffl)) . The actua l value of y 
depe nds on the galaxy populati on considered and fo r submillimeter galaxies is estimated to be /i ma x ~ 10 — 30 (|Perrotta et al 



2002 



3). However, as discussed in Paciga et al. 1 200&t l. there are still large uncertainties on the sizes of SM Gs, and there is one 



SMG lensed by the Abell 2218 cluster, whose magnification is estimated to be [i ~ 45 ( Kneib et al. 20o4 ) 



5.1.3 Schechter-hke distribution 



Since the effect of lensing depends on the local slope of the differential counts, we expect much larger magnifications if the 
intrinsic distribut ion is sufficientl y steep. For illustrative purposes now we assume that the intrinsic counts distribution is of 
a Schechter type ( Schechter 19761 ) 



*W _ (N\(s_y e _ s/S ^ (49) 



dS \S' J \S 

such that, after lensing, it becomes close to the measured BLAST distribution. Beyond the turnaround flux S' , this function 
becomes extremely steep and lensing is very effective. As an example we fix parameter values for this distribution as 7 = 2, 
N = 3.5 x 10 3 , S' — 3.5 mjy, and show the lensed distributions in Fig. (|11|) . Here the Schechter function fits intermediate 
points and decreases fast enough at faint fluxes so as not to exceed the CIB flux. The lensed distribution is consistent with 
the BLAST data given the error bars if these galaxies are at redshifts z s > 2. In this context, observed counts above ~ 40 mjy 
have all been highly magnified from intrinsic lower fluxes. Note the agreement with the brightest point can be improved at 
the expense of a larger disagreement with intermediate points by e.g. increasing the normalization factor N. 

In the absence of strong theoretical motivation for the intrinsic distribution, and given the relatively large error bars 
on the faint end, we did not attempt to fit any particular model to the observed counts. However, given the two extreme 
cases considered here, it can be seen that the data can be fit by a variety of models, with very different implications for the 
evolution of the luminosity function of these galaxies. 
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Figure 12. Blackbody intensity Iqmb (thick solid line) and SZ flux densities AS 1 ^ 2 for halos at zi 
10 14 2 to 10 1B h^M® (dashed and dot-dashed lines). 



0.5 and virial masses ranging from 



5.2 SZ Mass Estimation 

The CMB has specific intensity spectrum 7cmb of a blackbody with temperature Tomb ~ 2.725-K". In terms of x = 
hv/k B T c ,MB = ^/(56.78GHz) and I = 2h / c 2 (k B T C MB / h) 3 = 2.699 x 10 8 Jy sr" 1 we have 

7cmb = B v = ^{e h "^' r ^-iy 1 = ^ l . (50) 

As CMB photons cross the cluster hot gas and interact with its high temperature electrons via inverse Compton scattering, 
the outgoing photons gain energy, shifting their occupation number, temperature and intensity. 

The change in specific i ntensity can be c omputed by solving the full relativistic Boltzmann equation and expanding in 
powers of 8 e — (ksT e / 'm e c 2 ) ( Itoh et alj|l998h to obtain 

^=f{x,T R )g{x)y, (51) 
where T e is the electron temperature, g(x) = x 4 e x /(e x — l) 2 , the Compton i/-parameter is defined as 

y = a T (9 e / n e dl = [ n e (l)T e (l)dl , (52) 



and f(x,T e ) — fo(x) = x(e x + 1)/ (e x — 1) — 4, neglecting terms of higher order in e , which are negligible in the non-relativistic 
regime. Since I v = dS v /dQ, the change in flux AS„ through the cluster is computed by integrating over the cluster solid angle 
dfl = dA/d\ 

AI v dn = I f(x,T e )g(x)Y, (53) 



with the integrated Compton F-parameter given by 



Y = ydQ= -^4 / n e T e dldn = / n e T e dV , (54) 

J m e c 2 J d\ m e c 2 J 

and the cluster volume element is dV = dAdl = d\dQ,dl. One can model electron density n e and temperature T e profiles and 
compute Y from them. Since the total number of electrons in the cluster N e = J n e dV is proportional to the total cluster 
mass times the gas fraction /gas, one expects the scaling j n e T e dV oc Af gas T gas . Further considerations for the case of a gas 
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in virial equilibrium produce the expected scaling relation 



We employ fits from hydrodynamic simulations of Nagai (2006) for this relation given by 



d A (z) J V 1014/l " 1M o, 

Therefore, given the cluster mass M and redshift zi we can estimate its integrated Y parameter through the scaling relation 
and the change in flux density follows from Eq. (1531) . 

In our results we consider two mass definitions for the scaling relation. We use parameter values A14 = 2.5, a m = 5/3, 
which are appropriate for the mass M v j r defined within the virial radius. These values are roughly intermediate between 
adiabatic simulations and those with cooling and star formation processes. Similarly, we take A14 = 7.0, a m = 5/3, for the 
mass M2500 defined within the radius r2soo where the halo overdensity is 2500 times the critical density. 

In Fig. [12] we show the blackbody spectrum 7cmb as well as the SZ distortion for halos at zj = 0.5 and with different 
virial masses. Here we employed the non-relativistic limit f(x,T e ) = fo(x) for illustrative purposes. 



5.2.1 Contamination of the SZ signal 

We now estimate the effect of SMGs and their magnification on the SZ signal of clusters. SMGs have typical spectra that 
peak at ~ 100/xm, but they still have significant flux around ~ 1mm, where the SZ effect becomes important. 

Since the SZ effect is measured by subtracting the flux of the CMB photons in the line-of-sight of the halo from the 
average flux of background CMB photons, adding an average flux S„ from background galaxies in principle does not affect the 
signal, since the background is removed by the subtraction. This remains true even with lensing since the average flux does 
not change (Eg. I44|) . However, local flux fluctuations a eal (S v ) in the background flux due to Poisson noise in galaxy counts 
cannot be removed. 

The true SZ flux A5^ z is then modified by a sal (S^). These fluctuations are further increased by lensing magnification to 
°-obs(SV) so that the 68% range of the observed SZ flux becomes 

(ASf )obs = ASf ±o%j,(S v ) (57) 

In the absence of lensing effects the flux noise is simply o-^ E (S v ) = a gal (S'^) = \J S 2 V where the fluctuation is within the 
cluster virial radius 

= AQ vil J dS„Sl^- . (58) 

In Fig. 1131 we show the SZ flux for clusters of various masses at zi — 0.5 and the contributions of background SMGs. Here 
the average flux and its noise are computed with the BLAST distribution at A = 500 /im (v — 600 GHz) and extrapolated 
to other frequencies using the SED prescription. The average background flux SV from these galaxies within the virial radius 
of the cluster is rather significant: it is larger than the SZ signal. However, as a "background" it can be subtracted. Without 
lensing, the noise a s ^(S v ) in the flux is substantially smaller than the background, but non-negligible in comparison to the 
SZ signal. 

For clusters of large mass, even though the fluctuations from background galaxies increase due to the larger cluster area, 
these fluctuations become a smaller fraction of the total SZ signal, and relatively less important. Less massive clusters however 
are highly affected by the fluctuations of SMGs. For clusters of mass 10 1 ,6 /i _1 M0, expected to be above the minimum mass 
detected by upcoming surveys, even the unlensed contamination is comparable to the SZ increment. The right panels of 
Fig. [13] show the effect of using a different aperture in the definition of mass and SMG contamination. Note the relative level 
of contamination remains about the same without lensing. 

In both panels, magnification due to the cluster is rather significant, especially in the case of SZ fluxes through smaller 
radii (r25oo), since those are confined to cluster cores of high magnification. The flux noise is magnified by factors of ~ 2 — 7. 
Magnification due to large-scale structure in the line of sight is much smaller - it enhances the noise by only a few percent. 

The left panel of Fig. 1141 sho ws results for clusters of virial mass 1O 14 ' 6 /i _1 M0, but with the Lagache model for the flux 



from SMGs (jLagache et al.ll2004r ). As mentioned earlier, the Lagache model agrees well with the BLAST data at A = 500/im, 
so the dotted and dashed lines are nearly identical to those of the middle left panel of Fig. 1131 For the Lagache model, the count 
distributions are available at wavelengths within the SZ range. The triangle and square dots shown here are computed directly 
from the Lagache model at A = 500, 850, 1380, 2097/xm {u = 600, 350, 220, 140GHz). Notice that the SED extrapolation from 
A = 500/^m (dotted and dashed lines) fits well the model points. 

The right panel of Fig. [14] shows the results for the BLAST data, but for a cluster redshift of zi = 1.0. The qualitative 
results again do not change. We have verified this for other masses and redshifts as well. 



© 0000 RAS, MNRAS 000, 000-000 



18 Lima, Jain and Devlin 





0.4 




0.2 







1 1 


-0.2 


>\ 


U. 1 


1-3 




' ' 


0.05 


U 

'> 





U 




V 


-0.05 


GO 






0.04 




0.02 









-0.02 



1 1 1 1 1 1 1 


1 1/ 1 1 

/ X 
s / 






1 


1 1 1 1 


M„ lr =10 15 M e /h: 
1 i i i 1 i i i 





1 1 1 1 1 1 1 1 1 


/ 

/ 


1 1 ,.-| 1 yV 1 | 1 1 1 - 


- i 


iii 


M vir =10 14S M e /h; 
1 1 1 1 1 1 1 1 1 1 



_ 1 


11,1 
3 / 

»/ 


i ,) • i ^ i i 


| i i i^ 

- Unlensed: 
■ Lensed 


- 1 


iii 


M vir = 

ill 


10 14 2 M Q /h£ 
1 i i i 




200 400 600 800 

v [ GHz ] 



200 400 600 800 

v [ GHz ] 



Figure 13. SZ flux for clusters of different masses at zi = 0.5 (solid thick line) and the corresponding average background flux S v and 
noise cr sa ^(Su) within the cluster radius from a population of SMGs at z B = 2.0 (dashed lines). (Left): Fluxes measured within the v irial 
radius r v j r . (Right): Fluxes measured within r25oo- Radii and masses are converted with the prescription of Hu & Kravtsov (2003) for 
a NFW profile. For clusters of larger masses, the fluctuations from background galaxies increase, but become a smaller fraction of the 
total SZ signal. Lensing magnification (dotted lines) further enhances the flux noise, especially within smaller apertures. 




Figure 14. SZ flux for clusters with a virial mass of 10 14 - 6 Ji _1 M o , as in the middle left panel of Fig. 1131 (Left): Instead of using the 
BLAST data, we employ the Lagache model for estimating contamination by SMGs. The dashed and dotted lines are still computed 
extrapolating from the v = 600GHz value with the SED prescription. Triangle and square dots denote average fluxes and Poisson 
fluctuations at exact wavelengths A = 500, 850, 1380, 2097/^m, and show that the extrapolation agrees well with the model. (Right): Same 
as middle left panel of Fig. 1131 but for cluster redshift Z; = 1.0. Both A5^ z and the contribution of SMGs decrease for higher 2;, but 
the contamination level remains roughly the same. 
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6 CONCLUSION 



We have implemented a halo model calculation of magnification effects of galaxy and cluster halos. We explored the reliability 
of our model and compared it with published results from simulations. The one-halo term computed in the model describes 
well the magnification probability at large magnifications, but und ersestimates it at regio ns of de-magnification (fj, < 1). 
The effects of compensating halos and addition of a two-halo term ( Coorav fc Sheth 2002h . as well as the effects of cluster 
substructure can further improve our semi-analytical calculation. 

We applied magnification probability distr ibutions from nume rical simulations (Hilbert et al. 2007) to the counts of 
distant galaxies. Our results differ from those of Paciga et al. ( 20081 ) : our expression differs from theirs by two magnification 
factors, leading to a significantly smaller lensing contribution than their estimate. For the counts of high redshift SMGs 
recently observed, we find that lensing leads to enhancements of at least 60 percent at the bright (steep) end of the counts . 
In particular cluster fields with large magnification cross- sect ion, this enhancement can be significantly larger. 

Galaxy cluster surveys that utilize the SZ effect are susceptible to at least three kinds of contamination: by radio galaxies 
associated with the cluster, by high redshift SMGs that lie behind the cluster, and by projection effects due to large scale 
structure. The latter may result for instance fr om additional group s and clusters below the survey detection threshold which 
end up being projected on observed clusters IjHallman et al.ll2007T ). Here we estimated the contribution of SMGs over the 
frequency range targeted by SZ surveys. We used the counts from the BLAST survey, extrapolated to lower frequencies, and 
included the effects of lensing magnification from simulations and semi-analytical estimates. 

We find that for clusters with masses between 10 14 and 10 15 solar masses, there is a significant contamination from SMGs 
at frequencies higher than the SZ null at 220 GHz. This contribution arises from the irreducible variance of the Poisson 
distribution of the SMGs. It is comparable to the SZ excess at 350 GHz at the low mass end, and is about half the SZ excess 
for a mass of M v ir = 4 x 10 14 solar masses (it can be either an excess or a deficit for a particular cluster). With arcminute 
sized beams, SMGs are confusion limited, so it is not possible to isolate and remove their contribution. These result s are 
consistent with earlier studies of the (unlensed) contribution estimated from SCUBA sources (jWhite fc Maium dai 2004). We 
also find that lensing from the cluster cores significantly enhances this irreducible noise, especially for SZ fluxes within small 
appertures containing the cluster critical curves. 

Our results imply that SZ surveys must model the contamination noise of SMGs at the SZ null and for channels at higher 
frequencies. The value of the signal measured in such channels is significantly contaminated, even for the highest mass clusters 
due to lensing magnification. At frequencies below the SZ null, the SMG contribution is at the few percent level for mass 
ranges of interest. So it is less likely to be a problem in cluster detection in the SZ decrement regime, but must be included for 
cosmological measurements, as even a few percent bias in the inferred cluster mass can affect derived cosmological parameters 
due to the steepness of the mass function (especially at the high mass end). 
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APPENDIX A: LENSING BY ELLIPSOIDAL HALOS 

The effects of lensing on source images are introduced in Section 3.1. It is useful here to recall basic results for the axially 
symmetric lens case before considering the elliptical case. In general the deflection angle vector a — (01,02) relating source 
and image angles can be obtained from the convergence field k(9) by 

ai ) = ij ( fe' K{ e')I^, (Ai) 

where 9 = (6i, 62)- If the density profile has axial symmetry, so do all statistical quantities derived from it. In particular the 
lensing potential is independent of the position angle with respect to the lens center f(9i, 62) = <p{9). Choosing the lens center 
as the origin of the coordinate system, the deflection vector points towards the lens center with magnitude 

a{9) = ±J = , (A2) 

where m{6) = M{< 9)/nDfE crit with M(< 9) being the mass within 9. 

Given an axially symmetric lensing potential <p(8), we obtain the elliptical generalization with major axis along the 92 
direction by replacing 9 2 — 9\ + 9\ by 



92 

1! L 02 



+ 9 2 {l-e v ). (A3) 



(l-e„ 

Using the derivatives 89 )89\ and 89/892, the Cartesian components of the deflection angle are given by 
8 V gi s 

ai = W^JT^jS^' 

- = H-y-^m, (A 4) 

and their derivatives are 

dai a{9) Bl f a{9) da(9)\ 

897 ~ (l~e^)9 + {l-e v ) 2 9 2 \ T + ~dW ) ' 
802 a(fl)(l-e„) e' 2 2 {l-e v ) 2 f a{0) d&(9) \ 

89 2 9 9 2 V 9 d9 ) ' 

8 ai 8a2_9 1 9 2 { a(9) da{9)\ 



892 89i 9 2 V 9 d9 

Here and below tildes denote spherically symmetric quantities, e.g. a(9) = dtp/d9. First, let us check the axially symmetric 
case by setting e v = 0, in which case 9 2 = 9 2 — 9\ + 9\ and 

ddi a{8) OiBj ( &{9) da{9)\ 

Wi - + ^ 2 -{-- r + ^ r ) ■ (A6) 

In this case, the convergence and shear simplify to 

l(a(9) da(9)\ 
2 \-^ + —) ' 

fl2 



-A(_a(e) da(9) \ . _ 6162 f a{6) da(0) 



71 29 2 V 8 + d9 ) ' 72 9 2 \ 9 + d9 



d9 

Notice that the average value of k within 9 is 

"® = "4 f WR{9')d9' = ^B-, (A8) 
tt9 2 J 9 
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and therefore it is easy to verify the fact that 

™ - K^-^H- 

Using the following equalities 

50) _ . . 

in the expressions for general ellipticity Eqs. (|A5|) we obtain 

dai 2 | fli 

901 (1-e^) l7l (l-e v ) 2 2 ' 

^- = (/t+|7|)(l-e y )-2|7| ^ 2 , 

Therefore the convergence is given by 

k = A(e v )K(0) + B(e v ,6>i,0 2 )|7(0)|: (All) 
where 

AM ~ 2(1 -e v ) ' B (<V.fc.fe)-^)- 2(1 _ e ^ )2 - 
and the shear components are given by 

7i = c(e v )k(e) + D(e^,e 1 ,e 2 )\^(e)\, (A12) 

where 

CM = ^n 1 "^ 2 . ^,6> 1 ,fl 2 ) = C(e y )- 6,? - 6,2(1 -^ )4 



2(l-e v ) ' ~^>— ' ~™ 2 (l- e „) 2 



and 



72 = -^|7(5)|. (A13) 

Therefore, we can compute k and 71 and 72, at positions 61, 82 from the values of k and \ j\ at 8 from the spherical case. 
The calculation of j —y | and p, then follows. 



APPENDIX B: HALO MASS CONVERSIONS 



Simulations of the SZ effect usually predict its flux decrement /increment for different cluster mass definitions. Observationally 
it may be desirable to define SZ halo masses at small apertures (large overdensities), since those can provide fluxes with larger 
signal-to-noise ratios. Moreover one expect smaller radii to enclose less flux from radio sources and SMGs (though less SZ 
signal as well). However, as we show, smaller apert ures also correspond to incr eased flux fluctuations from lensing. We follow 
the mass conversion recipe outlined in White ( 2001 '! and Hu fc Kravtsov ( 2003h . which we briefly review here for overdensities 
defined in terms of the critical density. 

Converting different ha lo mass definitions re quires assuming a prescription for the halo profile. Throughout, we have been 
assuming the NFW profile (jNavarro et al.lll997l ) 



p(r 



M h 



(r/r 3 )(l + r/r s ) 2 ' 

For this profile, the mass Mh within radius rn enclosing an overdensity Ah with respect to the critical density p c 



4-jrr 



-A h pc = 4np s rlf(r B /r h ) . 



where the first equality is a definition, the second results from the NFW profile and 
f(x) = z 3 [ln(l + a:" 1 ) - (1 + x)- 1 ] . 



(Bl) 



(B2) 



(B3) 
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For the usually employed virial mass M v \ r the general equation above becomes 



4-Kr 



— A c pc = 47rp s rv ir /(l/c) 



Mvir 

where c = r v i T /r s is the halo concentration. From these equations we can relate M% to M v i r as 
M h A h 1 



M vi: 



A c (cx)3 ' 



= r 



i / A h 



/(1/c) 



^(/) = 



We employ the fitting formula for inversion of x (/) = / _1 (|Hu fc Kravtsovll2003h 

1/2 

+ 2/, 



«i/ 2p + ( 7 



(B4) 

(B5) 
(B6) 

(B7) 



where p = a 2 + a 3 In / + a 4 (ln f) 2 and (ai,...,a 4 ) = (0.5116,-0.4283,-3.13 x 10 _3 ,-3.52 x 10 -5 ). This fit is accurate to 
better than 1% for the range of virial masses and concentrations used here. 
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